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Abstract 

We analyze how temporal variability in local demography and dispersal combine to affect the rate of spread 
of an invading species. Our model combines state-structured local demography (specified by an integral or 
matrix projection model) with general dispersal distributions that may depend on the state of the individual 
or its parent. It allows very general patterns of stationary temporal variation in both local demography and 
in the frequency and distribution of dispersal distances. We show that expressions for the asymptotic spread 
rate and its sensitivity to parameters, which have been derived previously for less general models, continue 
to hold. Using these results we show that random temporal variability in dispersal can accelerate population 
spread. Demographic variability can further accelerate spread if it is positively correlated with dispersal 
variability, for example if high-fecundity years are also years in which juveniles tend to settle further away 
from their parents. A simple model for the growth and spread of patches of an invasive plant (perennial 
pepperweed, Lepidium latifolium) illustrates these effects and shows that they can have substantial impacts 
on the predicted speed of an invasion wave. Temporal variability in dispersal has gotten very little attention 
in both the theoretical and empirical literatures on invasive species spread. Our results suggest that this 
needs to change. 

Keywords: Invasions, invasive species, spatial population dynamics, stochastic demography, fluctuating 
environment, integral projection model, perennial pepperweed 



1. Introduction 

Invasive organisms are "altering the world's natural communities and their ecological character at an 



unprecedented rate" (Mack et al. 2000 p. 706), and often have substantial impacts on ecosystem structure 
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and function (SimberlofF 2011). Troublesome invasives in our current home states include zebra mussels, 
emerald ash borer, giant hogweed, avocado thrips, and smooth cordgrass (Spartina) . 

Spatial models of population spread have a potentially important role in evaluating and designing strate- 



gies for preventing or slowing the spread of invasive species (e.g., Taylor and Hastings (20041; Grevstad 



( |2005[ ) ; [Jongejans et al.| ( |2008[ ) ; [Bogich et al] ( |2008[ ) ; |Epanchin-Niell and Hastings] ( |2010[ ) ) . For well-studied 
species, simulation models allow detailed demographic modeling and accurate representations of landscape 
structure (e.g., Higgins et al. (2000); Jongejans et al. (2008); Andrew and Ustin ( |2010 ); Minor and Gardner 



(2011)). But simple "strategic" models (such as deterministic or stochastic matrix models) have often been 
useful for identifying the life stages or demographic processes that are the best targets of opportunity for 



management efforts to preserve a native species or control an invasive (e.g.. Shea and Kelly (1998); Heppell 



et al. (2000); Morris and Doak (2002); Lande et al. (20031; Shea and Kelly (2004); Shea et al. (2010)). 



Sensitivity analysis of the long-term population growth rate A, or of the long-term population spread rate 
c* , have often been a key tool in these applications. 



In a seminal paper Kot et al. ( 1996 1 showed how integrodifference equations could be used to model 
realistic patterns of organism redistribution (e.g., long-tailed distributions rather than the Gaussian spread 
that results from classical reaction-diffusion models), and gave a simple expression for the asymptotic rate 
of population spread. They found that long-tailed dispersal distributions can give faster rates of spread than 
a Gaussian distribution with the same mean square displacement. This has been proposed as a resolution 



of "Reid's Paradox", the rapid northward advance of tree species after the last glacial retreat (Clark 19981. 



The analysis by Kot et al. ( 1996 ) was quickly extended to include temporal variability in local population 



growth (Neubert et al. 2000), demographic stochasticity (Lewis 2000), discrete stage structure (Neubert 



and Caswell 2000), and two-dimensional spread (Lewis et al. 2006). Two recent extensions have been 



models with continuous population structure (Jongejans et al. 2011) and models that combine discrete 



stage structure and demographic variability (Schreiber and Ryan, 2011 Caswell et al. 2011 1. Caswell et al. 



(2011) also provide formulas for sensitivity analysis of the long-term population spread rate for periodic or 
stochastic environmental variation. 

Here we take two additional steps. The first is the natural step of combining demographic variability 
with continuous population structure, using the formalism of integral projection models. The second is to 
analyze the effects of an ecologically important aspect of population spread that has received surprisingly 
little attention in previous work: temporal variability in dispersal. Empirical evidence is very limited, but 
suggests that the frequency and range of long-distance dispersal can vary greatly from one year to the next 



(Andrew and Ustin 2010). The general formulas for population spread rate and its sensitivity in models 



with temporal variability (Neubert et al. 2000 Schreiber and Ryan 2011 Caswell et al. 2011) allow for 
temporal variability in both local demography and dispersal, but previous to this paper there has not been 
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(to our knowledge) a mathematical analysis of how dispersal variability can affect the rate of population 



spread. However, a recent numerical study by Seo and Lutscher (20111 did examine periodic fluctuations 



in dispersal rates for populations with sedentary and randomly diffusing individuals. They found that 
fluctuations in dispersal rates could increase or decrease rates of spatial spread depending on temporal 
correlations between dispersal and demography. To better understand these and other interactions between 
demography and dispersal on population spread, we analyze three forms of variability, separately and in 
combination: temporal variation in local demographic parameters (e.g. survivorship and fertility), a single 
mode of dispersal whose parameters (e.g., mean dispersal distance) vary over time, and multiple modes of 
dispersal (e.g., local wind dispersal and long-range animal dispersal) whose frequencies vary over time. 
Our analysis reveals that the effects of dispersal variability can be very different from those of demographic 



variability. A classical (Lewontin and Cohen 1969) and very general result (Tuljapurkar 1990 Ellner and 



Rees 2009 ) is that temporally uncorrelated demographic variability reduces population growth and spread 



rates (Lewis 2000 Clark et al. 2001 Schreiber and Ryan 2011 Caswell et al. , 2011). In contrast, we find 



that temporally uncorrelated dispersal variability can increase the rate of population spread. Moreover, when 
dispersal is variable rather than constant, the effect of demographic variability can be reversed: demographic 
variability that by itself would decrease population growth and spread rate can instead increase those rates, 
if it is correlated with dispersal variability. These general results are all derived by perturbation analysis for 
small fluctuations, but we also provide a simple geometric explanation for the effect of dispersal variability. 
We then use an empirically-based model for the spread of an invasive plant (perennial pepperweed) to 
show that our results continue to hold at very high levels of variability and that dispersal variability and 
dispersal-demography covariance can have appreciable effects on population spread rate. 



2. Model and assumptions 

We consider a continuously structured population in which the state z of an individual (e.g. size or 



age) lies in a compact set of all possible individual states Z (Ellner and Rees 2006). These individuals 



disperse along a one dimensional transect of their environment (however, rates of spread in a two dimensional 
region can be computed by "marginalizing" a two-dimensional dispersal kernel along the direction of interest 



(Lewis et al. , 2006)). Consequently, the location x of an individual can be identified with a point on the line 



X = (—00,00). Let nt{x,z) denote the population density at location x, state z, and time t. In the absence 
of density dependence (which we will consider in section HI), the most general form of the model is 



nt+i{x, z) 



Kt{x, z, xo, zo)nt(a;o, zo)dxodzo 



(1) 



where x is location, z is individual state, nt{x, z) is the population distribution in space and state at time 
Kt is the kernel for year t, and the integral runs over the spatial domain X = (—00, 00) and the (compact) 
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set of possible individual states Z. The kernel Kt{x, z, xq, zq) represents the rate at which individuals in state 
zq and location xg at time t produce individuals in state z and location x at time t. It includes changes in 
individual state, changes in location, and production of new offspring which may vary in state and location. 
We will often write 7if+i — KtUt as a shorthand for equation ([T]), and similarly for other kernels. 



Consistent with prior studies ( Kot et al. 1996 Neubert and Caswell 2000 Neubert et al. 2000 Jongejans 



et al. 2011 Caswell et al. 2011 Schreiber and Ryan 20111, we assume spatial homogeneity. In particu- 



lar, state transition rates are the same at all locations, so our model incorporates temporal variability in 
environmental conditions but not spatiotemportal variability, and movement probability is a function of the 
distance between the starting and ending locations. That is. 



Kt{x,z,xo,zii) = Kt{x- a;o,z,zo), with Kt{v,z,zo) = Kt{-v,z,zo). 



(2) 



Within that constraint, however, the dispersal pattern can depend on individual state in any way, in principle. 
Any constraints on movement dictated by the species' life history is reflected in the structure of the kernel. 
For example, if new offspring undergo natal dispersal (e.g., seeds or larvae) but then settle for the rest of 
their life (e.g., trees, corals), the kernel has the form 



Kt{v, z, zo) = SQ{v)Pt{z, Zq) + kd.t{v)Ft{z, Zo) 



(3) 



where 5q is the Dirac delta-function (a unit mass at u = 0), and P are the fecundity and survival/growth 
kernels respectively, and kd^t is the juvenile dispersal kernel that describes the displacements of offspring 
from their parent. Without loss of generality we assume that kd^t is a probability distribution, i.e., that any 
offspring mortality prior to establishment is absorved into F . 

To ensure that invasion speeds are well-defined, we need several additional assumptions. First, we assume 
that dispersal events have exponentially bounded tails. More precisely, we assume that the transformed 
kernels 

Hs,t{z,ZQ)= j Kt{v,z,zo)e'''"dv (4) 

are finite with probability 1 for all s in some interval (— si,si); the interval is symmetric because of our 
spatial homogeneity assumption. In the case of juvenile dispersal, equation ([s]), the transformed kernels are 



H,^t^Pt + Mt{s)Ft 



(5) 



where Mt is the moment-generating function of kd^t- In equations Q and ([s]), s characterizes the shape of 
the invasion wave, and the kernels Hg t determine the spread rate for an invasion wave where total population 



density decreases exponentially at rate s as a function of distance from the population center (Appendix A I 



Without an exponentially bounded tail, the rate of spatial spread may constantly accelerate, so there is no 



asymptotic invasion speed (Kot et al. 1996). 
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Second, we assume that for all s in (— si,si), the temporal sequence of transformed kernels 



{. . . Hg^-ij Hsfl, Hsi, . . .} arc stationary, ergodic and satisfy the assumptions of EUner and Rees (2007) 
for stochastic integral projection models. Stationarity means that the pattern of temporal variability in the 
environment doesn't change over time: each year is (randomly) different from the previous, but there are no 
long-term trends in the statistical properties of the variability, such as an increase in rainfall or temperature. 
The other main assumptions are that the kernels must be continuous, bounded, and in an ergodic set. The 
ergodic set property is analogous to the requirement that a matrix projection model be irreducible and 
aperiodic. It means that for some m > 0, i?s,m ■ • • ^s,2^^s.i is positive at all (z, zq) with probability one, 
where the kernel product is defined as 



H2Hi{z,zn) = J H2{z,y)Hi{y,zo)dy 



representing the effect of Hi acting on the population followed by H2. For the juvenile dispersal kernel (|3|, 
the transformed kernels will satisfy this assumption if the base kernels do and the temporal variability in 
Mt{s) is bounded (—00 < mi < Mt{s) < 777,2 < +00 with probability 1). 

Under these assumptions, the dominant Lyapunov exponent 7(3) of the random sequence of transformed 
kernels exists, i.e. there is a number 7(5) > such that 

lim l\og\\H,,tHs,t-i---Hs.i\\=j{s) (6) 



with probability one. Here, |j-ffs,t|| denotes any operator norm of H. 



s,t- 



3. Population grov^rth and spatial spread in the linear model 

We begin by examining population growth and spatial spread in the absence of density-dependent feed- 
backs. This analysis provides key information about how the population grows and spreads when rare. In 
particular, when the population experiences negative density-dependence, it determines the rate of spatial 
spread at the leading edge of the invasion and hence the asymptotic invasion speed (see section |4]) . 

To understand asymptotic population growth without density-dependent feedbacks, we integrate both 
sides of ([1]) over space and arrive at the overall state distribution ht{z) defined as 



nt{z) = J nt{x,z)dx 

So long as it is finite, this overall state distribution changes according to the demographic kernel Kt{z, zq) = 
Ho,t{z,zo) Kt{v,z,zo)dv, i.e. 

nt+i{z) ^ Kt{z,zo)ht{zo)dzn. (7) 
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The total population summed over all locations is therefore described by a nonspatial stochastic model. So 
starting from a finite initial population, the total abundance exhibits asymptotically exponential growth 
or decay. Whether the population grows or shrinks is determined by the dominant Lyapunov exponent 
7(0) of the the random kernels Kt- In the ecological literature the dominant Lyapunov exponent is usually 
called the "stochastic growth rate" logAs where A5 is analogous to the dominant eigenvalue of a matrix 
projection model. Under our assumptions, if the population starts from a finite initial distribution, then 
limt_s.oo 7 In / / nt{x, z)dxdz — log As with probability one. So if A5 > 1 (equivalently, if 7(0) > 0) the 
population grows exponentially, and if A5 < 1 it shrinks exponentially, in the long run. If the environmental 
variability is not too strongly autocorrelated, then the probability distribution of total population size is 
approximately lognormal, and the mean and variance of log-transformed total population size are both 



linear functions of time (EUner and Rees 2007) 



When As > 1, the total population size increases and one would expect the total area occupied by the 
population to increase. In an infinite-population model such as ours, "occupied" has to be defined by a 
threshold density nc, so that location x is regarded as occupied at time t if the total population (or total 
biomass, or some other measure of total population size) at location x at time t is at least ric- When the 
initial spatial and state distribution of the population is of the form n{x, z) = 710(2)6^''^ with s > 0, we show 
in [Appendix A| that the asymptotic wave speed is 

( \ ^^^^ fSi\ 

c{s) = — (8) 

where 7(5) is the dominant Lyapunov exponent of the transformed kernels i?s,i- But unlike these initial 
distributions, real invasions begin with populations whose densities are bounded and restricted to a finite 
spatial interval. For these realistic invasions, we show in [Appendix A| that the asymptotic wave speed is 
bounded above by 

c*= mill c(s). (9) 

Therefore, c* is the only biologically relevant wave speed. 

The relevance of c* hinges on the conjecture that an initially localized invasion will asymptotically spread 
as a wave, and therefore at rate c*. Simulations for unstructured and matrix population models universally 
support this conjecture for dispersal kernels with tails that decrease exponentially or faster, so that a finite 
moment generating function exists. However, this has been proved rigorously only for unstructured models 
under some additional assumptions about the shape of the dispersal kernel (see MoUison 1991). 



Our result generalizes prior results on wave speeds in fluctuating environments. Neubert et al. (2000) 



found the same formula for unstructured populations living in a temporally uncorrelated environment. More 



recently, Caswell et al. (2011) and Schreiber and Ryan (2011) independently derived this formula for matrix 



models (i.e. the kernel K is supported by a finite number of states) in random environments. Schreiber 



and Ryan (20111 also showed that the average spread rate up to time t is normally distributed with mean 



c* and a variance that decreases inversely proportional to t. The results of Ellner and Rees (2007) imply 



that the same conclusion applies to the integral projection models considered here, when the environmental 
variability is not too strongly autocorrelated. 

4. Density-dependent models 

To account for density-dependent effects on dispersal and demography, we allow the transition kernel for 
individuals at location xq to depend on the local population density n{xo). Then the most general form of 
the model is 

nt+i{x,z) ^ J Kt{nt{xo);x,z,xo,Zo)nt{xo,zo)dxodzo (10) 

where the kernel Kt{nt{xo)', x, z, xq, zq) depends on the local density at xq- Despite this additional complexity, 
there is compelling evidence that the asymptotic invasion speed is determined by the linearization at n = 0. 
This broad principle is called the "linearization conjecture" and is expected to hold provided per-capita 



survivorship and reproduction are greatest at low densities (MoUison 1991). This conjecture has been 



mathematically verified for structured populations exhibiting compensating density dependence and living 



in constant or periodic environments (Lui 1989 Weinberger 2002). Moreover, there is growing numerical 



support for this conjecture for structured populations exhibiting overcompensating density-dependence and 



living in constant or fluctuating environments ( Neubert and Caswell , 2000 Neubert et al. 2000 Schreiber 



and Ryan 2011 Caswell et al. 2011). 



For the general nonlinear model ( 10 ) considered here, the linearization conjecture applies provided two as- 



sumptions are met. First, population growth and dispersal is greatest at low densities i.e. Kt(0; x, z, xq, zq) > 
Kt{nt{xo); x, z, a;o, zq). Second, the kernel for the unoccupied habitat Kt{0; x, z, xq, zq) satisfies the assump- 
tions discussed in the previous section i.e. spatial homogeneity i^t(0; Xq + v, z, xq, zq) — ^t(0; z, zq), and 
the existence, stationarity and ergodicity of the transformed kernels Hs^t{z,ZQ) ~ J Kt{0;v, z, Z())e~''^^dv. 
The first assumption is violated for populations exhibiting positive density dependence at low densities (an 
AUee effect). For such populations, whether the population spreads or not can depend on the initial size 



and spatial distribution of the population (Kot et al. 1996). When spatial spread occurs, these populations 



grow and spread faster when they achieve higher densities. Consequently, the linearization conjecture can 
underestimate asymptotic invasion speeds for populations exhibiting an AUee effect. 

5. Variable dispersal and spread rate 

A consistent finding in integrodifference models of population spread has been that temporally uncorre- 



lated demographic variability slows the rate of spread (Neubert et al. , 2000 Caswell et al. 2011 Schreiber 



and Ryan 2011). The root of this phenomenon is the smaU-fluctuations approximation for the stoehastic 



growth rate 7 — logAg, see, for example, (B.5). To leading order, any kind of temporally uncorrelated 



variability about average demographic rates causes the population growth rate to decrease in proportion to 



the interannual variance (Lewontin and Cohen 1969 Tuljapurkar 1990 Ellner and Rees 2007). 



In contrast, as we will now show, temporal variability in the range of dispersal can increase the rate 
of spread. A heuristic explanation for this difference is that Jensen's inequality acts in opposite directions 
on demographic and dispersal variation. For demographic variation, population growth depends (roughly) 
on the geometric mean of annual population growth rates, or equivalently on the mean of log-transformed 
annual growth rates. The logarithm function has negative second derivative, so variability decreases the 
mean. Dispersal variation affects the spread rate through M{s), the moment generating function of the 



dispersal distribution. M{s) has a positive second derivative (see Appendix C), so variation in dispersal 



can increase the mean of the kernels Hs_t that determine the spread rate. In many situations the increase in 
mean dominates the negative impact of the increased variability due to dispersal variation, so spread rate 
increases. 

This is the most technical section of the paper. Because our main conclusions have been stated in the 
preceding paragraph and the Abstract, readers who wish can avoid the technical details by skipping from 
here to section [6l 



5.1. Sensitivity formulas for invasion speed 

Our results are based on analyzing the effects of small perturbations to the kernel Kg^t- A small (order 
e) perturbation of K causes a perturbation to the transformed kernels from Hs.t to Hs_t + £ Csj for some 
generally random kernel C. To approximate the effects of this perturbation using Taylor series, we need 
derivatives with respect to e at e = 0. 

Two general results which apply to any perturbation are 



dc* 1 

oe s* 



1 c'"7 
s de" 



{Vs\t+l,Cs',tWs',t) 



(11) 



(12) 



where E denotes the expectation of a random variable, (/, 5) = //(x) • g{x)dx denotes the inner product 
of two integrable functions, the invasion speed c* associated with the base kernel Hs^t occurs as s = s*, and 
Vs,t,Ws^t are the time-dependent stationary reproductive value and population structure sequences of the 
base kernel H^.t- Some important perturbations (such as adding variability in dispersal distances) affect the 
"noise" structure but have no mean effect, i.e., E[Cs^t] — 0. Provided that the Hg^t and Cs,t are independent 
and identically distributed, and Cs^t is independent of -ff^,*, the effect of these random perturbations on 
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invasion speed is of order 



— — = and = E 

oe oe^ s* 



(13) 



The effect of the perturbation on asymptotic wave speed is therefore times the right-hand side of (13 1 



Equation ( 13 ) also holds if Cg t is not independent of the unperturbed process but its mean, conditional on 



any event in the unperturbed process, is identically zero. 



The derivations of (111, (12 1 and (13) are in Appendix B Methods for numerically evaluating the 



expectations in (12) and (13) are given by Ellner and Rees (2009). 



5.2. Single dispersal mode 

As a simple illustration of how variable dispersal can accelerate spread, consider the juvenile dispersal 
model ([3]) when the juvenile dispersal kernel has constant shape but varying mean dispersal distance. That 
is, kd^t{v) — {1/Lt)ki{v/Lt) where fci is a dispersal kernel satisfying our assumptions in which the mean 
absolute parent-offspring displacement is one unit of distance. The scaling by Lt guarantees that k^^t is a 
probability distribution, and Lt is then the mean absolute parent-offspring displacement in year t. Let Mi(s) 
be the moment-generating function of ki; then 



Hs,t = Pt + M,iLts)Ft. 



(14) 



For small fluctuations in Lt, we let Lt = L + a^zt where E[zt] — 0, Var[zt] = E[z^] = 1 and ctl ^ 1. We 
assume that the zt are independent and identically distributed, and independent of the demographic kernels 
Pt,Ft. Using the small- variance approximation 



Mi{Lts) = Mi{Ls + saLZt) « Mi(Ls) -t- saLM[{Ls)zt + 
the transformed kernel is 



Ft + M,{Ls)Ft +saLM[(Ls)ztFt + '-^^l^^i^z^Ft. 



(15) 



(16) 



The first two terms on the right-hand side of ( 16 ) are the constant-dispersal transformed kernel which we 



denote and the last two terms are the perturbation due to variable dispersal. 



The last term on the right-hand side of ( 16 ) has positive mean. Using the sensitivity formula ( 12 ) and 



the independence of zt, the increase in the mean of the transformed kernel due to this term increases the 



invasion speed by (see Appendix D for details 



.s*alM['{Ls*) 



E 



{Vs*,t,FtWs',t) 
{vs',t,H^,tWs',t) 



(17) 
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where c(s*) — c* is the hwasion speed for the unperturbed kernel -ff^ and Ws»,t, Ws*,t are the time-dependent 
stationary reproductive value and state distribution for H^, j. Opposing this effect is a decrease in the 
invasion speed caused by the kernel variance. The dominant contribution to the kernel variance comes from 



the Zt term on the right-hand side of (16 1, which has zero mean. Using the perturbation formula (13) with 



(Tl as the small parameter e, this term decreases the invasion speed by 



s*al{M[iLs*)y 



E 



{Vs*.t,FtWs*.t) 



(18) 



Variable dispersal will increase the spread rate whenever (17) is greater than (18). Because ff^j > Mi{Ls)Ft, 
we have 



< 



1 



{Vs,t,FtWs,t) 



(19) 



Therefore (17 1 is greater than (18) if Af"Mi > (M()^. We show in Appendix C that this condition always 
holds, at all points where the moment generating function is finite. Thus, small variation in the range of 
dispersal (added to variability in survival, growth and fecundity) always increases the asymptotic rate of 
spread in the juvenile dispersal model. 



In Appendix F we generalize this result to state-dependent dispersal, i.e. a dispersal kernels kd^t{v) 



{\/ Lt {z, zo))ki{v/ Lt{z, zq); z, zq) so that movement distributions can depend on both the parent and offspring 
states. For example, taller plants would be expected to have longer average seed dispersal distances. If the 
temporal variance in Lt is small, and in any one year Lt is either above-average for all (z, zq) or below-average 
for all (z,zo), then variance in Lt increases the asymptotic spread rate. 

While this simple example illustrates that temporal variation in mean dispersal distances can increase 
the asymptotic invasion speed, not all forms of dispersal variation increase invasion speeds. Moreover, there 
can be subtle, nonlinear interactions between forms of dispersal variation that individually have opposing 
effects on invasion speeds. To illustrate these opposing trends and nonlinear interactions, we examine next 
populations with multiple dispersal modes. 

5.3. Multiple dispersal modes 

Organism frequently use different modes of dispersal. These modes may corresponds to different envi- 
ronmental currents (e.g., marine propagules at different heights in the water column experiencing different 
water currents), multiple vectors (e.g. plants dispersed by bird, mammal or insects), or individual differences 
in physiology or behavior. Let ki (v) be the dispersal kernel associated with the i-th mode of dispersal and pi 
the probability of an individual utilizing this mode of dispersal. Then the population level dispersal kernel 
equals k{v) = '^i^iPiki{v) and the associated moment generating function is M{s) — '^i^i Pi'mi{s) where 
rrii is the moment generating function for fcj. 
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We consider two forms of fluctuations on these multiple modes of dispersal: small fluctuations in the 
frequencies p{t) — {pi{t), . . . of these different modes of dispersal, and small fluctuations in the mean 

dispersal distance Li{t) for each of the modes of dispersal. For simplicity, we focus on the case where there 
are no fluctuations in local demography i.e. Ft = F and Pt = P for all t. The dispersal-demography kernel, 
in this case, equals Hs,t = P + ^^Pi{t)m.i{L,i{t)s)F . Define 

£i{t) ^ Li{t) - L^, Qi{t) = Pi{t) - pi, 

where Xj = E[_Li(i)] is the expected mean dispersal distance for dispersal mode i, pi is the expected frequency 
of dispersal mode i, and s* determines the invasion speed c* — c*{s*) for the unperturbed kerned i/^ = 
P + pimi{Lis)F . If we assume the fluctuations are small (of order e), a Taylor expansion to second order 
of the kernel Hg t at s — s* yields 

Hs'^t ~ P + ^p,(t)(m, + m,,i^,(t)+m,^2^»(i)')^' 

i 

« P + ^p^m.F + ^ {Q^(t)mi + p^m^,l£^{t)) F + ^(ft(t)mi,i^,(<) + p^m,^2^i{tf)F . 



H" 



The perturbation At has mean zero and variance of order e^. The perturbation Bt has non-zero mean of 
order and variance of order and, consequently, this variance is negligible. 

Let Aq, w, and w be the the dominant eigenvalue, the reproductive value and the stable state distribution 



for iJg. , respectively. Using the sensitivity formula (12) and the fact that {v,H^,w) = Aq, the correction 
term to the invasion speed c* due to the perturbation Bt is 



1 



s*Ao 



mv,Btw) 



(w, Fw) 
s*Xo 



(m,,iCov[ft(i), £,(*)] -Kp»m,,2Var[£,(t)]) . 



(20) 



Since the perturbation At has zero mean, the sensitivity formula ( 13 ) implies that the correction term to 
the invasion speed c* due to the perturbation At equals 



1 



jE[{v,Atwf] = - 



{v, Fw)' 



Var 



Yg^{t)mi + pim,^iii{t) 



We can expand the variance term as 



Var 



Var 



2Cov 



Qi {t)mi, ^ Pifn^^iii {t) 



(21) 



(22) 
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The correction terms (20) and (21) provide several insights into the separate and combined effects of 
variation in the frequencies of the dispersal models and the mean dispersal distances. Since the variance 



Var (3i(t)mi] in the frequencies of dispersal modes only show up in the second correction term (21 1, 
small fluctuations in these frequencies per-se decrease invasion speeds. All else being equal, negative cross- 
correlations between dispersal modes wth similar kernels increase invasion speeds (i.e., if one mode of long- 
range dispersal is rare when another is common), while positive cross-correlations between similar modes 
decreases invasion speeds. 



The variances Var [£i{t)] in mean dispersal distances appear in both correction terms (20) and (21 ). The 
combined effect of these terms on the invasion speed equals 



(w, Fw) 
s*Xo 



^Pj?Tii^2Var[£j(t)] 



(w, Fw) 
2Ao 



Var 



y^^Pimi^i£i{t) 



(23) 



We show in Appendix C that (23) is always positive. Therefore, the net effect of the variability in mean 



dispersal distances is to increase the invasion speed. Moreover, equations (22) and (23) implies that negative 
cross-correlations between the £i{t) result in a larger positive effect on the invasion speed c*, while positive 
cross-correlations ameliorate this effect. 

The net effect of fluctuations in the frequencies in the dispersal modes and their mean dispersal distances 



is determined by their separate effects plus the covariance terms in equations (20) and (21). The combined 
effect of these covariance terms on the invasion speed equals 



(u, Fw) 
s*Xo 

{v, Fw) 



^TOi,iCov[ft(t),^j(i)] 



{v, Fw) 



ft(0 - 



Ao 

V, Fw) 



Gov 



gt{t)mi,'^Pimi^il^{t) 



Ao 



-Pi 



(24) 



While (24) shows that generally the correlations between frequency {pi{t)) and range {£j(t)) of dispersal 
potentially have complex interactive effects on the invasion speed, this interplay simplifies significantly when 
the rrii are all equal e.g. all dispersal modes have the same expected mean dispersal distance and underlying 
movement kernel. In this case, the fact that J^jPji^) = 1 implies that J^j Qj{t) — ^-nd (24) simplifies to 

(u, Fw) 



s*Xo 



E"^.,iCov[P,(t),L,(i)] 



(25) 



Thus, faster rates of spatial spread occur if a given dispersal mode has a higher mean dispersal distance 
when it is more common. 

To illustrate some of the implications of these approximations, we consider the simple case of an unstruc- 
tured population exhibiting two modes of dispersal. The frequencies of these models are given by p{t) — pi (t) 



12 



Li(t) = 1 ±0.01 and LsO) = 2±0.02 



Li(t) = 1±0.9and L2(t) = 2±1.9 




Figure 1: Invasion speeds for an unstructured population with two modes of dispersal. Dispersal frequencies given by p(t) = pi {t) 
which takes on the values 0.5 it cr with equal probability, and P2{t) = 1 — p{t). The mean dispersal distances are given by 
Li{t) = Li ± Ti with equal probability. The correlation between Pi{i) and Li(t) is p, while the correlation between Pi(t) and 
Lj{t) with j 7^ i is -p. In A, n = 0.01 and T2 = 0.02. In B, n = 0.9 and T2 = 1.9. In both figures, P = 0.6, F = 1.5, Zi = 1, 
L2 = 2, and 'mi{s) are Laplace moment generating functions. 

which takes on the values 0.5 ± cr with equal probability, and P2{t) — P{^)- The mean dispersal distances 
are given by Li{t) = Li±Ti with equal probability. The correlation between Pi{t) and Li{t) is p, while the 
correlation between Pi{t) and Lj{t) with j ^ i is —p. Fig. illustrates, as asserted by approximation (21 1, 



that larger fluctuations in the frequencies of the dispersal models decreases the invasion speed provided 
that fluctuations in the mean dispersal distances are sufficiently small. In contrast, when these fluctuations 
are sufficiently large and positively correlated with the corresponding dispersal mode frequencies. Fig [Tj3 
illustrates, larger fluctuations in the dispersal modes can increase invasion speeds. 

5.4- Interactive effects of demographic and dispersal variation on spatial spread 

Temporally uncorrelated fluctuations in demographic rates (e.g. survivorship, fecundity), as discussed 
earlier, reduce the stochastic growth rate log A5 of a population and, as a consequence, reduce invasion speeds. 
However, when there are simultaneously fluctuations in demographic rates and mean dispersal distances, 
correlations between demography and dispersal may alter this conclusion. To investigate this possibility, 
consider the juvenile dispersal model with one mode of dispersal and assume that the fluctuations in the 
mean dispersal distance Lj, the fecundity kernel Ft, and the survivorship and growth kernel Pt are small. 
Define 

TO = Afi(Ls*), mi ^ M[{Ls*)s*, tos = M{'(Zs*)(s*)V2, 
it = Lt-L, TTt^Pt-P, (j)t^Ft-F, 
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where L = E[it] is the expected mean dispersal distance, F — E[i^f] is the expected fecundity kernel, 
P = E[Pt] is the expected survival/growth kernel, and s* determines the invasion speed c* = c*(s*) for the 
unperturbed kerned H'^ = P + Mi{Ls)F. If £t, ttj, and (j)t are small (i.e. of order e), a Taylor approximation 
of the moment generating function up to order yields 



H 



Pt+Ah{LtS*)Ft 

Pt + (m + mitt + ma^t ) Ft 



P + mF + TTt + rrKpt + mi 



At 



Bt 



The perturbation term Bt has non-zero mean of order and variance of order and, consequently, this 
variance is negligible. The perturbation term At has mean zero and variance of order e^. 

Let Ao, w, and w be the the dominant eigenvalue, the reproductive value and the stable state distribution 



for iJg. , respectively. Using the sensitivity formula (12) and the fact that {v,H^tw) = Ao, the correction 



term to the invasion speed c* due to the perturbation Bt is 



1 



1 



AqS' 



E[{v,Btw)] ^ — {miCoy[£t,{v,(j)tw)]+m2{v,Fw)Ya.T[it]). 



(26) 



Since the perturbation At has zero mean, the sensitivity formula ( 13 ) implies that the correction term to 
the invasion speed c* due to the perturbation At equals 



1 



1 



„ ^„-E[{v,Atw)^]^ ^ 
2s* A^ ' 2s* A^ 



Var [{v, TTtw) + m{v, (ptw) + miit{v, Fw)] 



(27) 



where we can rewrite the variance term as 

Var [(w, TTtw) + m{v, (ptw) + mi£t{v, Fw)'\ =Var [{v, TTtw) + m{v, (j)tw)] + Fw)'^Yai [it] 

+ 2mi{v,Fw) {Cov[{v,TTtw)Jt] + mCov [{v,(l)tw),£t]) 



(28) 



The correction terms (26 1 and (27) provide several insights into the separate and joint effects of variation 
in dispersal and demography on invasion speeds. The variance in the mean dispersal distances Var[£t] occurs 
in both correction terms (26) and (27|. However, the fact that Af"Afi > (M{)^ implies, just as argued 
in Section 5.2, that the net effect of dispersal variation per se is to increase invasion speed. Demographic 



variance occurs only in the second correction term (27). Consequently, the effect of variance in survival or 
fecundity per se is to decrease the invasion speed by 

1 



2A^s^ 



:Var [{v, TTtw) + m(f , (j^tw)] 



(29) 



This effect is largest if survival and fecundity are positively correlated, and smaller if their correlation is 
negative. 
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The covariance Cov[(w,7rti 



between survival and mean dispersal distance occurs only in the second 



correction term (27). Consequently, all else being equal, a positive dispersal-survival covariance decreases 



the invasion speed, while a negative correlation increases the invasion speed. The covariance Cov[(w, (ptw) , it] 



between fecundity and mean dispersal distance occurs in both correction terms (26 1 and (27 1. The net effect 
of this covariance from both correction terms is 



mi 



^^^)-^Cov[(.,<^,.),4](«,P.) 
Ao / s*Ao — ^ — ' 

>o 



(30) 



Demographic variability therefore makes two contributions: one negative (equation (29l) and the other 



positive when fluctuations in fecundity and dispersal are positively correlated (equation (30l). If the second 
of these is the larger in magnitude, the net effect of demographic variability is to increase invasion speed, 
a phenomena demonstrated numerically in our analysis of a patch-based model of pepperweed spread (see 
section 7 and Fig|3]). Equation (30) also implies that for annuals (i.e. P = 0), the covariance between 
dispersal and fecundity has no effect on the invasion speed. In fact this observation holds in general in 
our plant model with juvenile-only dispersal, not just for small fluctuations. With Hs^t = Mi{Lts)Ft, the 
Lyapunov exponent 7(5) is the sum of ElogAfi(Lts) and the dominant Lyapunov exponent for {Ft}. Hence, 
covariance between Lt and Ft has no effect on the invasion speed for annuals. 



6. Is dispersal variability important? 



Andrew and Ustin (2010) used remote-sensing data to estimate four interannual dispersal kernels for 



patches of perennial pepperweed, Lepidium latifolium, in the Sacramento-San Joaquin River delta in Cali- 
fornia, and developed a simulation model for pepperweed spread. Dispersal differed among habitats and was 
also highly variable from year to year. While most new recruits established close to their parents in all years, 
long-range dispersal (> 100 m) occurred in two interannual periods (« 0.2% and 2% of new recruits). As a 
result, the average dispersal distance (combining all sites and habitat types for each year) had a temporal 



coefficient of variation of approximately 40% (Andrew and Ustin 2010, Table 3). 

However, they noted that "invasion dynamics are little affected by temporal variation in dispersal dis- 
tances. This contradicts theoretical expectations from analytical models, which predict that temporal vari- 
ation in either reproductive rate or dispersal will reduce spread rates to the geometric mean of the rates 
observed under the component constant conditions (Neubert et al. 2000)." Specifically, simulated rates of 
spread with temporally variable dispersal were very similar to the rate when the dispersal kernel with the 
farthest long-distance dispersal was used in all years (their Figure 6). This appears to to contradict our 
analytic results above on the importance of dispersal variation, as well as previous theory. 



In this section, we explain why the findings of Andrew and Ustin (2010) are consistent with our results 



and with previous theory, and actually are predicted by our results. The explanation is a phenomenon that 
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Figure 2: Numerical calculations of wave speed for the unstructured perennial plant population described in the text, with 
parameter values Li = 1, F = 1.5, P = 0.6. (A) Asymptotic wave speed c* as a function of p2i the probability of a long-distance 
dispersal year, for L2 = 5, 10,20 (see legend of panel B). The curves show c* for 0.001 < P2 < 0.95; the solid circle is c* for 
P2 = 0, which does not depend on the value of L2. (B) c* / L2, the asymptotic wave speeds relative to the scale of long-distance 
dispersal. (C) s*L2, the asymptotic wave-shape parameters corresponding to c*. (D) The graph of ■y{s)/s for L2 = 5 and 
various values of p2- The minimum values of •y{s)/s, indicated by horizontal dashed lines, correspond to the asymptotic wave 
speed c*. The vertical dashed lines correspond to the positive vertical asymptotes of f{s)/s when P2 = and when p2 > 0. 



has been observed numerically in models with temporally constant dispersal: rates of spread are far more 
sensitive to the range of long-distance dispersal than to the fraction of long-distance dispersers (see Figs. 13 



and 14 in Neubert and Caswell (2000)). The corresponding phenomenon for temporally variable dispersal is 
that the rates of spread are insensitive to the fraction of years in which long-distance dispersal occurs, and 
are largely determined by the mean dispersal distance in those years. 

In Figure |2] we illustrate this phenomena for an unstructured population with juvenile dispersal. The 
transformed kernel is Mt{s)F + P where F is the per-capita fecundity, P is annual survival of adults, and 
Mt{s) is the moment generating function for the dispersal kernel in year t. The dispersal kernel is a Laplace 
(two-sided exponential) kernel with two possible distance parameters L2 > Li representing long- and short- 
range dispersal years. The temporal variability was uncorrelated, the distance parameter in each year being 
Li with probability pi and L2 with probability p2 = 1 — pi. Figure [2]'V shows that the invasion speed 
c* changes discontinuously at P2 = 0: a minuscule chance of long-range dispersal (in these computations, 
P2 = 0.001) is very different from no chance at all. But a hundred-fold further increase in the chance of 
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long-range dispersal (from 0.001 to 0.1) has relatively little effect, especially when L2 S> Li. Figures [2]B 
and[2p show that to a very good approximation, the wave speed c* and the corresponding value of s* that 
minimizes 7(s)/s are linearly proportional to L2 (so long as p2 > 0), but are only weakly dependent on p2- 

The geometric underpinnings of this discontinuity, as illustrated in[2|l), are threefold. First, Mt{s) is 
defined on the interval (0, 1) with probability one for p2 = 0. However, when there is a positive probability 
P2 > of long distance dispersal with L2 = 5, Mt{s) is defined with probability one only on the smaller 
interval (0,0.2). Second, when p2 — 0, the minimal value of 7(s)/s occurs at s* « 0.6 > 0.2, while with 
P2 > the minimizer of 7(s)/s must occur in (0, 0.2). Finally, 7(3) for p2 > is greater than 7(5) for p2 = 0. 
Collectively, these properties imply a discontinuous increase in the invasion speed at p2 = 0. 

Although we have used an unstructured model for simplicity and illustrative purposes, this phenomenon 
is quite general. For example, suppose that i/|j(z,zo) with i = 1,2 are general transformed kernels of 
the form Q such that HI ^ are well-defined with probability one on the maximal interval (— Si,Si) with 
si possibly infinite. Let Hs^t equal with probability pi and i with the complementary probability 
P2 = 1 ^ Pi- If li{s)/s is minimized at s\ > S2 and 72(5) > for < s < S2, then there is a jump 

discontinuity in the invasion speed at p2 — 0. Moreover, a lower bound for the size of this jump discontinuity 
is 7i(s2)/'S2 ~ / ^i- Since this argument relies on §2 being finite, these jump discontinuities do not 

occur for thin-tailed dispersal kernels such as a Gaussian dispersal kernel. 



7. Application: a patch-based model for pepperweed spread 



To illustrate our results and the role of variable dispersal we develop and study a simple model for the 
spread of perennial pepperweed. Pepperweed is a Eurasian crucifer that has spread widely in the western 
United States since it was introduced into the US in the early 20th Century. It is now invasive in wetlands and 



riparian zones throughout the western US (Leininger and Foin 2009), and an agricultural weed of concern to 



alfalfa and native hay growers (Blank et al. 2002 ). The model is based on population studies in the Cosumnes 



River Preserve (Hutchinson et al. 2007, Viers et al. 2008) and nearby sites in the Sacramento-San Joaquin 



River delta, California (Andrew and Ustin 2010). Pepperweed has been identified as a significant threat 



in the Preserve because it is highly invasive there and forms monospecific stands where native species are 
completely eliminated. Our model incorporates the main features of pepperweed's life history and dispersal, 
but given the limited information on some demographic processes we can only claim that the model is inspired 
by pepperweed. The main purpose of the model is to explore interactions between temporal variation in 
local demography and dispersal. 
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7.1. Pepperweed life history 



Rosettes develop in the spring and eventually bolt to produce multiple inflorescences 0.5-2m tall (Leininger 



and Foin 2009). Seed production can be quite high, with high germination and viability (roughly 3000 



seeds/inflorescence and 174,000 seeds/m^ within a patch in the Preserve, with 96.4% germination and over 
80% viability seven months after seed production (Leininger and Foin 2009| )). However seed production and 



viability were substantially lower in moister or more saline sites ( Leininger and Foin 2009 1 



The root system of established plants spreads laterally and produces new sprouts, creating a dense local 
patch ft'om which all other species are eventually excluded. In suitable habitat, rhizomes can spread 1-2 



meters within a year (Forman Orth et al. 2006), and a single plant can grow into a patch several meters in 



diameter within 2 years (Leininger and Foin 2009). Root fragments can also disperse and resprout to create 



new patches. Patches are difficult to eradicate mechanically because of the deep root system and resprouting 
of new plants from root fragments. Pepperweed also forms a long-term seed bank, which can re-establish a 



patch after all adult plants have been killed (Viers et al. 2008). 



Water - soil moisture and flooding - has been identified as the main environmental driver of demographic 
variability. The study sites are seasonally flooded each year. In the Cosumnes River Preserve, established 
patches tend to shrink during wet years (i.e., growing seasons following a wet spring), and grow during dry 



years (Hutchinson et al. 2007). Patches can disappear in a wet year, but reappear in a subsequent dry 



year. At the "Experimental Floodplain" site in the Preserve, 153 patches (out of 443 total) disappeared 
between 2004 and 2005 and 40 (out of 312 total) disappeared between 2005 and 2006, but 75% of those 



had reappeared by 2007 (Hutchinson et al. 2007). New patches appeared mainly during dry years, but it 



is believed that these were initiated by seeds or root fragments that had been moved by water flow in wet 



years (Hutchinson et al. 2007). However, the effect of rainfall and flooding is site-dependent. During the 
same time period, at Bouldin Island nearby in the Sacramento-San Joaquin river delta, establishment of 
new patches was greatest in wetter years, with the largest effect in sites expected to be most water-limited 



(Andrew and Ustin 20101 



7.2. Model 

Because pepperweed occurs mainly in dense monospeciflc stands, we model the growth and spread of 
patches rather than individual plants. This allows us to use a density-independent model for the growth 
and spread of patches. Even though most plants experience strong intraspecific competition from others in 
their patch (even at the leading edge of the invasion) , a density-independent model of patch dynamics can 
be constructed for the situation where patches are sparse and inter-patch competition is rare. This model 
becomes inaccurate once patches are so big and common that they start to collide. But according to the 
linearization hypothesis, the linear model that applies to sparse patches at the leading edge of the invasion 
correctly predicts the asymptotic spread rate. 
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We assume the basic "plant" model, equation ([s]), approximating patches as expanding circles. Estab- 
lished patches have been reported to increase in diameter by 1-3 m/yr (Andrew and Ustin 2010), which is 
comparable to the initial rate of expansion from a founding plant. We therefore assume that patches grow 
at mean rate g{t) meters/yr, where g{t) depends on environmental conditions in year t, but not on patch 
size. To allow variance in growth among patches, we assume a lognormal distribution of growth rates with 
variance fT^, assumed to be constant. The log-transformed radius of new patches was assumed to follow a 
beta(3,3) distribution shifted to have mean -3; because of the rapid patch growth, the initial size distribu- 
tion has essentially no effect on model predictions. Patches "die" by disappearing in a wet year and never 
reappearing, and we let d(t) denote the annual probability of patch "death" (in the Experimental Floodplain 
data, E[d] « 0.06). 

We assume a Laplace (bi-exponential) dispersal kernel with time- varying mean dispersal distance L(t). 
This was based on the summary data in Table 3 of Andrew and Ustin (2010), which showed that for each 
habitat-year combination the mean and standard deviation of dispersal distance are roughly equal, and the 
observed maximum distance is roughly what would be expected under an exponential distribution given the 
mean distance and the sample size. We assume that L(t) has a mean of 15m (roughly the average over 
all years and habitats in Table 3 of Andrew and Ustm] (20101). We model one-dimensional spread, e.g., 



expansion along a strip of riparian habitat. 

We assume that fecundity - the founding of new patches by propagules from an established patch - is 
proportional to patch radius squared (i.e. to patch area), with time- varying constant of proportionality f(t). 

Let n{p, X, t) be the density at location x in year t of pepperweed patches of size p, where p is the natural 
log of plant radius. The demographic kernels making up the model are then 

Ft{p,po) = f{t)e^P" 13^,3(3 + p) (31) 

Pt{p,po) = {l~d{t)mp;\og{e'">+g{t)),a'g) 

where </>(•; /i, cr^) is the Gaussian density with mean p and variance cr^. The nonspatial kernel describing 
total population growth is then Kt = Ft + Pt, and the transformed kernels that determine the asymptotic 



spread rate is Hg^t = (1 ^ L{t)'^)^^ Ft + Pt- Because the model (31) is formulated in terms of continuous 
growth kernels, the demography is described by just a few parameters, which is very helpful for analyzing 
the model's behavior. 

7.3. Pepperweed model results 

We adjusted the mean of the fecundity coefRcient f{t) to match the observed average "birth rate" between 
2002 and 2007 in the Experimental Floodplain site (0.3 new patches per established patch per year). This 
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Figure 3; Numerical calculations of wave speed for the perennial pepperweed model. Mean values of the model parameters were 
L = 15,3 = 0.5,(1= 0.06,/ = 0.04 [og was held constant). (A) Temporal variability consisted of equally likely wet and dry 
years characterized by the fractional range of variability V, meaning that the two possible values of a parameter 9 are d± V9 
(9 = L,^,d, and/or /). "Vary demography" (solid black curve); L is constant, all other parameters vary in parallel: rapid 
growth, high fecundity, and low mortality during wet years. "Vary dispersal" (solid gray curve): L varies (further dispersal in 
wet years), all other parameters are constant. "Vary both": all parameters vary; with correlation C = 1, wet years are good for 
both local demography and long-distance dispersal, and with correlation C = — 1 wet years are good for long-distance dispersal 
but bad for local demography. (B) Temporal variability with four year types, characterized by all possible combinations of high 
or low long-range dispersal (L = L it VL) and good or poor local demography (e.g., g = g± Vg) with V = 0.8. C is the degree 
of correlation between fluctuations in local demography and dispersal. Solid curve: fecundity parameter / is considered to be 
a part of local demography. Dashed curve: / is associated with dispersal, so / and L vary in parallel. 



indirect estimate is supported by the fact that the model then predicts a long-term population growth rates 
of \s = 1.20 — 1.24 depending on the level of stochasticity, which matches very well the the observed rate of 
increase in the total number of patches in the Experimental Floodplain site ((456/158)^/^ ~ 1.24). 

Figure |3]4. illustrates the effect of dispersal variability alone and in combination with demographic vari- 
ability. These calculations used a simple wet years/dry years pattern of environmental variation, with the 
fractional level of variation V running from to 90% of the maximum possible (e.g., V — 0.5 means that the 
two (equally likely) possible values of the fecundity parameter / are f±0.5f, and similarly for the patch mor- 
tality d, patch growth rate 7 and dispersal range parameter L). Variation in dispersal alone (only L varies) 
increases the asymptotic spread rate, as predicted by our perturbation analysis, while variation in local de- 
mography alone decreases the spread rate (L constant, all other parameter varying). The small fluctuations 



approximation (equation B.5) implies that small temporally uncorrelated local demographic variation always 
decreases the spread rate, regardless of how demographic parameters covary. For our pepperweed model this 
remains true even when the difference between wet and dry years is very large. 
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More interesting is the interaction between the two types of variation. If local demographic variation 
is positively correlated with dispersal variation, so that wet years are good for both local patch growth 
and long-distance dispersal, then the effect of demographic variation is reversed, as predicted by equation 



(30): demographic variation now increases the spread rate c*, rather than decreasing it. Conversely, if the 
correlation between local demography and dispersal is negative, local demographic variability decreases c*. 
In Fig. [3] A, patch fecundity was regarded as a component of local demography, so when C = — 1 the years 
when propagules go far also tend to be years with fewer propagules. If we instead regard fecundity as a 
component of dispersal (so that L and / vary in parallel), adding local demographic variation on top of 
dispersal-distance variation still can increase the spread rate, but the effect of correlation is reversed: spread 
rate is increased by the addition of demographic variation when C = — 1, and decreased when C = 1. Positive 
correlation between variation in local demography and dispersal generally increases the spread rate (Figure 
[3]B), regardless of whether / is associated with local demography (i.e., high in years when survival and 
growth are high) or with dispersal (i.e., high in years when L is high). Spread rate is higher when fecundity 
and dispersal distance vary in parallel. This conforms to the general principle that spread rate is largely 
determined by the very best years for long-range dispersal, so making those years even better increases the 
spread rate even though overall fecundity is unchanged. 

Figure |4] show contour plots of the elasticity surfaces for population growth rate As and wave speed c* in 
response to changes in the mean value of a size-specific fecundity F or survival/growth P, for a wet /dry years 
scenario with negative correlation between local demography and dispersal (wet years are good for long-range 
dispersal, but bad for survival and fecundity). For example, the value of e{c* ^ F) at sizes pf pi,pt^i p2 
is the fractional change in c* per fractional change in the rate at which size pi parents produce offspring 
in the size range [p2 , P2 + dp] ■ Numbers in the bottom- left corner of each panel are the total elasticity (the 
integral of the plotted surface). For both A5 and c*, the survival elasticities are higher, a typical pattern in 
long-lived organisms because an increase in adult survival leads to a large increase in lifetime reproductive 
success. Total elasticities are higher for c* than for A5. The total elasticity of As is necessarily near 1 
(in a deterministic model the sum of the P and F total elasticities for A is exactly 1, here it is slightly 
different because of stochasticity) . The higher total elasticity for c* indicates that a small change in the total 
population growth rate produces a larger change in the spread rate. 

The relative unimportance of fecundity is less severe for c* than for As (panels C and D). This is again not 
unexpected - because it is only new offspring that disperse in this model - but it depends on the correlation 
between dispersal and local demographic variation. When we assume instead that wet years are good for 
both local demography and long-range dispersal, the relative importance of fecundity and survival/growth 
is about the same for c* as it is for As (Figure [s]). 

We also computed the elasticities of As and c* to perturbations of the mean values of the underlying 
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Figure 4: Elasticities with respect to mean kernel entries for the perennial pepperweed model with wet and dry year-types. 
The parameter values were the same as in Figure|3] with fractional range of variation V = 0.7, and correlation C = —1 so that 
wet years are good for long-range dispersal but bad for local demography (including fecundity), relative to dry years. The four 
panels show the elasticity surfaces e{y,x) for the stochastic population growth rate A5 and the wave speed c* with respect to 
the mean values of F{y,x) and P(y,x). Elasticity values are indicated by the contour lines, and by shading from zero (white) 
to high (dark grey). Numebers in the bottom-left corner are the total elasticity, i.e., the integral of the plotted surface. 

model parameters /, g, d, dg and the range of variation V (Figure |6|. Again, c* is more sensitive than A5 to 
model perturbations, and the most important parameter by far is the mortality rate of established patches, 
d (the graph shows the elasticity to a proportional decrease in d, so the effect is positive). The parameter 
elasticities are mostly insensitive to the pattern of correlation between local demography and dispersal, with 
two exceptions: the elasticity to / was higher, and the elasticity to V was lower, in the case that wet years 
are good for dispersal (high L) but bad for growth, survival, and fecundity. Increased variability V always 
increases spread rate, but has very little effect on population growth rate. The positive effect of variance 
in patch growth reflects the nonlinear relationship between patch size and fecundity, so that the gain from 
above-average patch growth more than offsets the loss from below-average growth. 
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Figure 5: The same as in Figure |4] but with correlation C = 1 so that wet years are good for both long-distance dispersal and 
local demography, relative to dry years. 



8. Discussion 

Stochastic fluctuations in demographic rates, such as survivorship, growth, reproduction, and dispersal 
are pervasive in natural populations and can have profound impacts on local population growth and persis- 



tence (Tuljapurkar 1990 Boyce et al. 2006 Ellner and Rees 2007 Benai'm and Schreiber 2009). Recently 



the implications of these fluctuations for rates of range expansion have been examined for discretely struc- 
tured populations ( Caswell et al.[ 2011 Schreiber and Ryan 2011). Here we have extended these studies, 
by allowing continuous structure (e.g., size) within populations, and by analyzing how interannual fluctua- 
tions in dispersal distance can affect the rate of spatial spread. Our paper continues the ongoing effort of 



extending the seminal work of Kot et al. ( 1996 ), who introduced to ecology the use of discrete-time integrod- 



ifference models for population spread with general (non-Gaussian) dispersal kernels. Using the formalism 



of demographic integral projection models (Easterling et al. 20001 Ellner and Rees 2006), we combined 



general models of local demography incorporating continuous and discrete population structure with general 
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Figure 6: Elasticities of population growth rate A5 and spread rate c* with respect to the mean of underlying model parameters 
/, g and Og (indicated by the symbol "s" in the plot). The parameter values were the same as in Figure |3] with two year 
types and fractional range of variation V = 0.7. In Case 1 (black symbols) wet years were assumed to be good for long-range 
dispersal and poor for local demography including fecundity /(<). In Case 2 (darker gray symbols) wet years were assumed to 
be good for long range dispersal and fecundity, and bad for patch growth and survival. In Case 3, wet years were assumed to 
be good for all demographic processes. The dashed line is the 1:1 line. 



dispersal distributions (potentially depending on the state of the dispersing individual and its parent), and 
general patterns of stationary environmental variation. 

For these general models, we show that the size of a species range asymptotically increases linearly 
with time. We provide an explicit expression for this invasion speed and, thereby, extend earlier results 



for unstructured populations (Neubert et al. 20001 and discretely structured populations (Caswell et al. 



2011 Schreiber and Ryan 2011). We also derive sensitivity formulas for these invasion speeds to small 



perturbations in demography and dispersal and, thereby, extend work of Caswell et al. (2011 ) to populations 
with continuous as well as discrete structure. From the analytic perspective, these sensitivity formulas allow 
us to explicitly investigate how small temporal variation in dispersal distances and local demography influence 
invasion speeds. From an applied perspective, these sensitivity formulas provide methods to evaluate how 
uncertainty in parameter estimates yield uncertainty in estimates of invasion speeds, and to identify which 
parameters have the largest impact on invasion speeds. 

Using the sensitivity formulas, we show that stochastic fluctuations in mean dispersal distances increases 
invasion speeds. In sharp contrast, temporally uncorrelated fluctuations in vital rates associated with always 
local demography decrease invasion speeds. This latter effect follows almost directly from classical results 



in stochastic demography (Lewontin and Cohen, 1969 Tuljapurkar 1990 Ellner and Rees 2007) where the 
population long-term growth rate is reduced by one-half of the "net" demographic variance experienced 



by the population (see e.g. equation (8) in Ellner and Rees (20091, which is the IPM version of equation 



equation (14.68) in Caswell (2001), Tuljapurkar's small variance approximation for temporally uncorrelated 
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environments) . Roughly, this reduction stems from population growth rate being a concave function of local 
demographic rates. Hence, by Jensen's inequality variation in these rates decreases the population growth 
rate and, hence, the invasion speed. In contrast, for the models considered here, the invasion speed is, roughly, 
a convex function of mean dispersal rates. Hence, fluctuations in mean dispersal distances increase invasion 
speeds. Intuitively, these fluctuations generate occasional bursts of long distance dispersal events and it is 
these occasional long distance dispersal events that effectively determine the invasion speed. This is consistent 
with the observation that dispersal kernels with "fatter tails" (e.g. Laplacian versus Gaussian) yield faster 



invasion speeds (Kot et al. 1996). Extending our analysis to reaction-diffusion models may explain a similar 



phenomena observed by Seo and Lutscher (2011). For populations with a sedentary reproductive stage and 



a randomly diffusing non-reproductive stage, Seo and Lutscher (2011) found that periodic fluctuations in 
diffusion rates also increase invasion speeds. 

Not all forms of temporal variation in dispersal enhance invasion speeds. Populations may exhibit multiple 



modes of dispersal either due to multiple dispersal vectors (polychory) (Berg, 1983 Bullock et al. , 2006 



Nathan et al. 2008 ) , multiple behavioral types within the population ( Fryxell et al. 2008 Cote et al. 



2010), or different passive forms of transport (Hampe 2004). The mean dispersal distances traveled in these 
modes can be substantially different. For example, Jordano et al. (2007) found that 50% of small passerines 
dispersed seeds of St. Luce Cherry less than 51 meters, while 50% of mammals dispersed seeds more than 
495 meters. Our sensitivity analysis reveals that temporal fluctuations in the relative frequencies of these 
modes of dispersal, in of themselves, reduce invasion speeds. 

The negative impact of fluctuations in the frequencies of dispersal modes or fluctuations in local demogra- 
phy can be reversed when these fluctuations are positively correlated with fluctuations in the "appropriate" 
mean dispersal distances. For example, if an increase in the frequency of a particular dispersal mode corre- 
lates with longer dispersal distances traveled by individuals in this mode, then variation in the frequencies 
of this dispersal mode can increase invasion speeds. Similarly, when variation in the local demography of 
individuals in a particular state positively correlates with the mean dispersal distance traveled by these 
individuals (e.g. fecundity covarying with natal dispersal distances), this variation increases invasion speeds. 
Hence, unlike correlations between local demographic rates which can only dampen the negative effects of en- 



vironmental fluctuations on population growth (Tuljapurkar 1990) and invasion speed, correlations between 
dispersal and demography can reverse the negative effects of fluctuations in local demography on invasion 
speed. We hypothesize that correlation between dispersal and local demography will be the rule rather than 
the exception. For example, larger plants will tend to produce more seeds, and release more seeds at a 
greater height so that they will tend to disperse further. Conditions allowing above-average growth for a 
plant species are therefore likely to result in higher fecundity and greater long-distance dispersal. 



For periodically- forced reaction diffusion equations, Seo and Lutscher (20111 numerically demonstrated 
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similar results for populations where individuals are born into a mobile, non-reproductive stage and mature 
into a sedentary, reproductive stage. They found that periodic fluctuations in birth rates could increase 
invasion speeds provided that these fluctuations were positively correlated with fluctuating diffusion rates. 
Interestingly, when individuals continuously move back and forth between mobile stages and reproductive 



stages, Seo and Lutscher (20111 found that negative correlations, rather than positive correlations, between 
birth rates and diffusion rates increase invasion speeds. 



Our results relate to the seminal work of Berg (1983) (cf. [Nathan et al. (20081) who emphasized the 
importance of 'chance dispersal' for long distance dispersal. Berg argued that chance dispersal comes in 
two differing forms. The first form involves "an unusually favorable combination of the regular dispersal 
factors" , occasions when many things happen to go well for dispersal at the same time. We have shown 
that positive correlations of dispersal distance with local demography or with frequencies of dispersal modes 
are "unusually favorable combinations" of factors that promote population spread. According to Berg, the 
second form of 'chance dispersal' corresponds to "an unusual coincidence involving a dispersal factor not 
normally operating together with the taxon in question." In our context, this "unusual coincidence" would 
represent a rare mode of long-distance dispersal. Our analysis reveals that even infrequent occurrences of a 
model leading to especially far-distant dispersal can produce a discontinuous and potentially large increase in 
invasion speed. This discontinuity has also been observed in temporally homogeneous models where a small 



fraction of the population continuously uses a mode of long distance dispersal. For example, Neubert and 



Caswell (2000) estimated spatial spread rates for a neotropical plant Calathea ovandensis whose seeds are 



dispersed by four ant species Horvitz and Schemske (19861. They found that the invasion speed is mainly 
determined by one ant species which disperses only 7% of the seeds but disperses the seeds the longest 
distance. 

We illustrated the synergistic effects of positively correlated fluctuations in dispersal distances and fe- 
cundity with a size-structured patch model of perennial pepperweed, an invasive crucifer of the western 
United States. When temporal fluctuations are sufficiently large and positive correlated with fecundity (e.g. 
correlated with preciptation), these correlated fluctuations increased invasion speeds by approximately 60% 
(cf . Fig. [3]) . Our sensitivity analysis of the pepperweed model also suggests that mortality of established 
pepperweed patches presents by far the best target for management. Even in the scenario that maximizes 
pepperweed spread (large, positively correlated variability in dispersal and demography), reducing the annual 
survivorship of established patches to 60% would push the population into decline and retreat (A5 « 0.8). 
Perennial pepperweed can be controlled by herbicide treatments, which achieve nearly 100% patch eradica- 



tion (Hutchinson and Viers 2011|, but non-chemical means are preferable or necessary in some settings (e.g. 



organic agriculture or natural areas). The results from our model support the conclusion of Hutchinson and 



Viers (2011) that some non-chemical methods would also be effective enough, such as a mow, till and tarp 
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treatment that was found to eradicate 50% of patches completely, and reduce by 95% or more an additional 
25% of patches. 

Collectively, our results illustrate the importance of fluctuations in dispersal kernels for rates of population 
spread. However, currently there are very few empirical studies which have measured these fluctuations, let 
alone their correlation with local demography. While estimating dispersal kernels is notoriously difficult, our 
work highlights that understanding a single snap shot in time of dispersal kernels may not be adequate for 
estimating invasion speeds. In particular, ignoring temporal variation in dispersal rates can substantially 
underestimate rates of spatial spread. From a theoretical perspective, our analysis raises the issue of whether 
temporal variation in patterns of dispersal alters other aspects of ecological dynamics, such as persistence or 
coexistence in spatially heterogeneous environments. Dispersal variation may also be important for processes 
other than species invasions, such as the spread of alleles (either new natural mutations or alleles escaping 
from a GMO crop) and the spread of novel pathogens. Because patterns of spread in alleles and rapidly- 



evolving pathogens such as flu can often be reconstructed through genetic analyses (e.g., Bedford et al. 



(2010)), the data needed for parameterizing models and testing predictions about spread rates might be 



most readily available for these applications. 
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Appendices 

Appendix A. Calculating the wave speed 

Let 7(s) denote the dominant Lyapunov exponent for Hs,f To find possible traveling waves, we consider 
population states of the form 

nt{x,z)=ut{z)e-'''. (A.l) 
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Assuming that this form holds at time t, we have (using the spatial invariance and symmetry of the kernel) 



nt+i{x,z)= II Kt{x — Xi, z, Zi)u{zi,t)e~^^'-dxidzi 

Kt{x — xi,z, zi)u{zi,t)e~'''^^^~'^^dxidzi 
Kt{v, z, zi)u{zi,t)e~^'"dxidzi 
Hs,t{z, zi)u{zi,t)dzi. 



(A.2) 



Therefore, if equation (A.l) holds initially, then by induction it holds for all time with Ut+i = Hg^tUt- The 



limit to the area occupied by the population at time t, Xt{s), is defined by the property that 

where w(z) is a non-negative weighting function (e.g., the biomass of a statc-z individual). Then 

1 



(s) = - log {Ut , w) - log 71c 



(A.3) 



and therefore 



Ms) , 7(s) 



as t — > oo 



with probability 1. In addition, if the environment process is uniformly mixing, then t^^ log (ui, w) has an 
asymptotic normal distribution with mean 7(3) and variance a2/t for some a. The average rate of spread up 
to finite time t, given by Xt{s)/t, is therefore asymptotically Normal with mean 7(5) and variance decreasing 
in proportion to . 



A standard monotonicity argument (see Kot et al. (1996) or Appendix A in Schreiber and Ryan (2011)) 



shows that any model solution with the initial population limited to a finite spatial domain cannot spread 
faster than the slowest such wave speed, given by equation (|9]). 



Appendix B. Sensitivity analysis of the spread rate 

Because c* is defined implicitly, for sensitivity analysis of c* we need to use the following general result. 
Suppose F{6) — min/(s, 0) where s, 9 are real and and / is a smooth real-valued function. Let s*{9) denote 

s 

the value of s at which f{-,9) is minimized. Then j^{s*{0),d) ~ and 



dF d 



86 de- 

Applying this to equation ([9]), we have 



df ds* df_ 



dc* 
1)6 



1 97(5*) 



(B.l) 



(B.2) 
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We can take to be a kernel entry Kt{v,y,x), an entry in one of the component kernels Ft, Pt, k^^t, or 
an underlying parameter in one of the kernels (e.g., the slope parameter in a regression model for survival 
probability as a function of size). A small perturbation to any of these has the effect of perturbing the 
transformed kernel Hs t to i/^.t + for e <C 1. 



There are two basic formulas for the resulting change in 7 (Ellner and Rees (2009); note that Ellner and 



Rees (2009) refer to 7 as logAs). The first, which applies to an arbitrary perturbation kernel Ct, is 

{Vs,t+l,CtWs^t) 



de 



= E 



(B.3) 



where Ws.t,w;s,t are time-dependent stationary reproductive value and population structure sequences of the 



base kernel Hg^t- The second concerns the situation where the right-hand size of (B.3) is zero, because the 
perturbation kernel has zero mean and is independent of the Hs^t- In that case, the change in 7 is O(e^) and 
to leading order 



7(s,e) = 7(s,0) 



-Var 



{Vs,t+l,Ct'Ws,t) 



= 7(s,0) 



-E 



{Vs^t+l,CtWs,t) 



(B.4) 



Applying (B.4 1 to perturbations consisting of small fluctuations around a constant mean kernel, we get the 

Var {vs,Hs,tWs) 



small-variance approximation, 



7 « log Ao 



2X1 



(B.5) 



where Aq is the dominant eigenvalue of the mean kernel E[_ffs^t] and Vs,Ws are the dominant left and right 
eigenf unctions of the mean kernel normalized so that {vs , Ws 



1. 



Ellner and Rees (2009) explain how 



these formulas can be used to compute the various sensitivities and elasticities that have been proposed for 
stochastic matrix and integral models, and the sensitivity to perturbations of the vital rate functions that 
are used to construct the kernel. 



Some of our results depend on the fact that (B.4 1 also holds under weaker assumptions than those 



stated above. In particular, the calculations in Ellner and Rees (2009) shows that (B.4) holds so long as 



the conditional mean of any perturbation kernel Ct given everything else (the full sequence of unperturbed 



kernels and of all other perturbation kernels) is identically zero. For example, (B.4) applies to perturbation 
kernels ZtFt where the zt are iid random scalars with zero mean and independent of the unperturbed kernels, 
even though this perturbation is not independent of the unperturbed kernel process. 



Appendix C. Moment generating functions of dispersal kernels 

We collect here some elementary but useful properties related to the moment generating functions of 
dispersal kernels. Let ^{dx) denote a probability distribution for spatially homogeneous dispersal with 
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displacement x, meaning that /i is a measure on the real hne and /i(^) = A). Let X denote a random 
variable with distribution /i and M{s) = E[e'*"^] its moment generating function. We assume that M{s) is 
finite on some (necessarily symmetric) open interval (— si,si). 

1. M(s) is infinitely differentiable on (— si,si) and its n*^ derivative is given by 

M(")(s) ^EiX^e"-^]. (C.l) 

Proof: This seems to be widely known but we have not found a proof in print. For s — this is a 
standard result about the moment generating function. The symmetry of fi implies that M{s) = M{—s) 
so it suffices to consider s > 0. We proceed by induction on n for an arbitrary fixed s e (0, si). For 
n = 1, M'{s) is the limit as e — > of 



E 



e 



(C.2) 



The integrand in (C.2 1 converges pointwise to Xe^ , as needed. To prove convergence of the integral 
we consider X > and X < separately, and find finitely /i-integrable upper bounds on the absolute 
value of the integrand, so that the Dominated Convergence Theorem applies. By Taylor's Theorem, 
the integrand in (C.2 1 is pointwise equal to Xe" for some s* (depending on X) between s and s + e. 
For X > 0, pick cti £ (s,si). Then for for aU e sufficiently smaU and X > 0, Xe"'^ < Xe"^^ = 
Xe~^^ e^'^^^^^^ . Choose 5 so that ai + 5 < si. For X > 0,Xe~^^ is positive and no larger than 



l/{6e). So for X > the integrand in (C.2 1 is bounded above by a constant multiple of e(°'i+'')^, which 
is integrable because cri + 5 < Si. For X < 0, the integrand is bounded above in absolute value by 
l/(s*e) which is uniformly bounded for e small. 

Assuming the result for n, the difference quotient for M("+i) is E [(X"e("+'^)^ - /e] . The 

construction of upper bounds for the integrand is nearly identical to the case n = 1, using the fact that 
^n+i^-Sx J > 0, a; > has a finite maximum. 
2. Unless Var[X] — 0, Af'^"'(s) is positive for n even, and has the sign of s for n odd. Proof: For n even 
this follows immediately from ( |C.l ). For n odd, because X and —X are identically distributed, we 



have 2Af(")(s) = E[X"e"^ + (-X)"e""^] = E[X"(e'*^ - e""^)], and the integrand has the sign of s 
pointwise except at X = 0. 
3. 7(5) is nondecreasing on (0, si). Proof: The kernel Kt{v, z, zq) can be factored as Kt{z, ZQ)iit{dv\z, zq) 
where fj,t{dv\z, zq) is the probability measure for state-z individuals at time t + 1 produced by state- zg 
individuals at time t. The last result, applied to the moment generating functions of the /it, implies 
that the transformed kernels H^^t are monotonically increasing in s on (0, si), therefore 7(5) must be 



everywhere nondecreasing. If the variability is temporally independent, then equation (B.3) implies 
that 7(5) is strictly increasing. 
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4. M'^'^^M > (Af(i))2 unless Yslt[X] = 0. Proof: using equation ( |C.1[ ), this is the Cauchy- Schwartz 
inequahty appUed to the functions c*^/^, Xe'^^/^, with strict inequahty because the two functions are 
not proportional to each other. 



5. Equation (231 is always positive. Proof: We will show that 



2Y,P^m^{L^s*) f ^ ) (s* )2 Var (f)] /2 - ^^^^ Var 



> (C.3) 



with a strict inequality if Var[£i(t)] > for some i. To simplify the notation, define 

Using this notation and noting that Gav[li{t),lj{t)] < yjYax[li (i)] Var[£j (t)] , we get that the left hand 



side of (C.3) is bounded below by 



(u, Fw) 
Ao 



E' 



(C.4) 



Using the fact that Aq = {v,H^,w) > {v,^^pimi{LiS*)Fw) ~ '}2iO'i{v,Fw), we get that (C.4) is 



bounded below by 



(C.5) 



Because m'l(LiS*)mi{LiS*) > {m[{LiS*)'^ as proved above we have that Oibi > Ci with strict inequality 
when Var[£i(i)] > 0. This observation plus the Cauchy-Schwartz inequality imply that 

E «^ E - (e ' ^ (e ' - ' > 

with a strict inequality if Var[^i(t)] > for some i. 



Appendix D. Two-parameter perturbations 



To understand the effect of the random perturbation ( 15 ) on the invasion speed, we initially consider the 



two parameter perturbation Hg.t + SiCg^t + £2^s,t where Hg^t — P + Mi{Ls)F, Cs.t = s{Lt — L)M[{Ls)F, 
and Dg^t — s^iLt — L)'^Mi{Ls)F/2. By Taylor's theorem, we have the second order approximation 



7(£i, £2, s) « 7(0, 0, s) + ^(0, 0, s)ei + ^(0, 0, s)e2 + ^-g^{0, 0, ^H + ^Q^iO, 0, s)e2 + ^^^(0, 0, 3)6^62 
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By sensitivity formula (12), we have 



9£2^ ' ' 2Ao(.) 
where (t£ — Var[_Lf] and Vs, Ws the corresponding left and right eigenfunctions of Hg scaled so that {vs 



1. By sensitivity formula (13), we have 
dj 



s^al{Mi{Ls))^{v,,Fwsr 



^^^(0,0,.)=Oand^(0,0,.)- 
Setting £i — e, e2 — e'^ , and only considering terms up to order e^, we get 



2Ao(s) 



7(0, 0,s) 



2Ao(s) 



M{'(Zs) 



2Ao(s)2 

(Af{(£s))^(i;,,fw, 
Ao(s) 



Appendix E. Sensitivity calculations for pepperweed 

We explain here how the general perturbation formulas were used to do the calculations for the pepper- 
weed model. The elasticity of As to the mean of a size-specific fecundity F{y,x) is by definition 



e{\s,F{y,x)) = 



F{y,x) dXs _F{y,x) dXc 



As dF{y,x) As dK{y,x) 
the last equality holding because K = P + F. From Table 3 in |Ellner and Rees (2009) we have 



ssiy,x) = AsE 



We therefore have 



In parallel. 



dK{y,x) 
eiXs,F{y,x))^Fiy,x)E 



vt+i{y)wtix) 
{vt+i,KtWt) 

vt+iiy)wt{x) 



vt+i{y)wt{x) 
{vt+i,KtWt) 



eiXs,P{y,x))^Piy,x)E 
The elasticity of c* to the mean of a size-specific fecundity F{y,x) is by definition 



(E.l) 



(E.2) 



(E.3) 



' dF{y,xy 

Because Hg.t = Mt{s)Ft+Pt, a unit perturbation in Ft{y, x) increases Hg^tiu, x) by Mt{s). The corresponding 



perturbation kernel Ct in (12) is therefore Mt{s)6y^x where Sy^^ is an approximate delta-function centered at 
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(jj,x), and the sensitivity at (y,x) is the hmiting value as the support of Sy^x shrinks to a point (EUner and 



Rees (2006)). We therefore have 



dc* 



dF{y,x) 



-E 



(E.4) 



Because c* = j{s*)/s*, we get 



e(c*,F(y,x)) = ^^^E 
7(s*) 



{Vs*,t+l,Hs'^tWs*,t) 



(E.5) 



For a perturbation to size-specific survival the factor Mt{s) is absent, so we get 



7(s*) 



{Vs*,t+l,Hs'^tWs 



(E.6) 



The expectations in the formula above can be computed for all (y, x) from one long simulation of the 



model using the methods described by Ellner and Rees (2009). For an IPM implemented using midpoint 
rule, these are the same as the standard methods for stochastic matrix models. For all such calculations we 
used 10000-year simulations, omitting "burn-in" periods of 250 years at the beginning and the end to allow 
convergence of the v and w vectors to their stationary distributions. 

For sensitivities to an underlying parameter such as the mean patch growth rate, the simplest approach is 
brute force: perturb the parameter by ±e, recompute wave speeds (reducing Monte Carlo error by using the 
same sequence of year types for all values of the parameter), and estimate the derivative by finite difference. 



Appendix F. State-dependent dispersal 

We consider here the plant model ([3| with dispersal dependent on the states of the parents and offspring, 
and small temporal variation in mean dispersal distance. Specifically, we assume 

kd.t{v, z, zo) = — rk{v/Lt{z, zo); z, zq) (F.l) 
Lt{z,zo) 

with {Lt~L) uniformly small and constant in sign as a function of (z, zq) for each t. The temporal variability 
in Lf is assumed to be iid and independent of the demographic variation. Note that the form of the dispersal 
kernel can also depend on parent and offspring states. Our goal is to show that the effect of the (small) 
variance in Lt is always to increase 7 and therefore to increase wave speed. 

We need the following general result. Let /, g, h, w be bounded nonnegative functions on a finite measure 
space {X, /i) such that fg > h? with strict inequality on a set of positive measure in the support of w. Then 
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letting ( , ) denote the L2 inner product on (X, fi) 



M = \ j V fwy/gwd^i 



(F.2) 



fwd^ij gwdfi 



The small- variance approximation to the mgf of the dispersal kernel (F.l ) is 



M{Lts) = M{Ls + {Lt ~ L)s) w M{Ls) + sM[{Ls){Lt - L) + —M"{Ls){Lt - Lf 



(F.3) 



with all terms in (F.3 1 being functions of (z,zo). Let o denote element-by-element multiplication, i.e. 



(/ ° 9){^j -^o) = f{z, ZQ)g{z, zq). The expansion of the transformed kernels Hg^t corresponding to (F.3) is 



Hs,t ~Pt + M{Ls) o Ft + sM'{Ls) o {Lt ^ L) o Ft + -M"{Ls) o {Lt - If o Ff (F.4) 



Because of the independence of Lt, the M' term in (F.4) has identically zero mean conditional on the 



unperturbed process. Its effect on 7(5) is therefore to leading order (omitting the s-dependence in u, w and 
M and its derivatives) 



-E 



{vt+i,{{Lt-L)oM'oFt)wt) 
{vt+i,H^wtf 



2n 



(F.5) 



where — Pt + M{Ls) o Ft is the unperturbed kernel. The M" term has nonzero mean, so its leading 
order effect on 7 is 



-E 



{vt+i,{{Lt^LfoM" oFt)wt) 
{vt+uH^wt) 



(F.6) 



We aim to show that (F.6) is larger than (F.5 1. We do this by showing that the integrand in (F.6) is with 



probability 1 larger than the integrand in (F.5), i.e.. 



{vt+i, {{Lt - Lf o M" o Ft)wt){vt+i,H^tWt) > {vt+i, {{Lt - L) o M' o Ft)wtY 



(F.7) 



Neither side of the inequality is affected if we replace Lt — L with its absolute value, so suffices to consider 
Lt{z, zq) — L > (recall that Lj — L is assumed to have constant sign in any given year). Because iJ" > MoFt 
and Vt+i,Wt are both everywhere positive, it suffices to show that 



{vt+i, {{Lt ~ If o M" o Ft)wt){vt+i, {M o Ft)wt) > {vt+i, {{Lt - L) o M' o Ft)wtf 



(F. 



Writing out the inner products as integrals, equation ( |F.8[ ) follows from (F.2) with / — {Lt — Lf o M" ,g 

M, h = {Lt - L) o M' and w(z, zq) = vt+i{z)Ft{z, ZQ)wt{zQ). 
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